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Abstract 

This work is devoted to the global stability theory of solutions for a general isothermal 
model of capillary fluids derived by C. Rohde in [17j , which can be used as a phase 
transition model. 

This chapter is structured in the following way: first of all inspired by the result by P.-L. 
Lions in [T3] on the Navier-Stokes compressible system we will show the global stability 
of weak solution for our system with isentropic pressure and next with general pressure. 
Next we will consider perturbations close to a stable equilibrium as in the case of strong 
solution. 

1 Introduction 

1.1 Presentation of the model 

The correct mathematical description of liquid-vapor phase interfaces and their dynam- 
ical behavior in compressible fluid flow has a long history. We are concerned with com- 
pressible fluids endowed with internal capillarity. One of the first model which takes into 
consideration the variation of density on the interface between two phases, originates from 
the XlXth century work by Van der Waals and Korteweg [13] . It was actually derived in 
his modern form in the 1980s using the second gradient theory, see for instance |12t I18j. 
Korteweg suggests a modification of the Navier-Stokes system to account additionally for 
phase transition phenomena in introducing a term of capillarity. He assumed that the 
thickness of the interfaces was not null as in the sharp interface approach. This is called 
the diffuse interface approach. 

Korteweg-type models are based on an extended version of nonequilibrium thermodynam- 
ics, which assumes that the energy of the fluid not only depends on standard variables 
but on the gradient of the density. In terms of the free energy, this principle takes the 
form of a generalized Gibbs relation, see [18] . 

In the present chapter, we follow a new approach introduced by Coquel, Rohde and their 
collaborators in [3]. They remark that the local diffuse interface approach requires more 
regular solution than in the original sharp interface approach. Indeed the interfaces are 
assumed of non zero thickness, so that the density vary continuously between the two 
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interfaces, whereas in the sharp interface models, the interfaces represent zone of dis- 
continuity for the density. Coquel, Rohde and their collaborators present an alternative 
model with a capillarity term which avoids spatial derivatives. The model reads: 



(NSK) 



with: 



dtp + div(pii) = 

d t (pu) + div(/ni (g> u) - pAu - (A + /u)Vdivu + V(P(p)) = KpVD[p] 
(p t =o,ut=o) = (j>o,uo) 



p > and A + 2p > 

where p denotes the density of the fluid and u 6 M. N the velocity, p and A represent the 
viscosity coefficients, At is a coefficient of capillarity, P is a general pressure function. We 
are particulary interested by Van der Waals type pressure: 



P : (0,6) -► (0,+oo) 

RT*p 2 
P(p) = ap 



where a, b, R, T* are positive constants, R being the specific gas constant. For fixed 
values a, b we choose the constant reference temperature T* so small as P to be monotone 
decreasing in some non-empty interval. 
Further we impose the conditions: 

u(t,x) — > 0, p(t,x) — >0 as |x| — > -)-cx), (1-1) 

In last section, we consider also more general situations: monotone pressure, other bound- 
ary conditions than boundary condition ll.il at infinity instead of ll.lt 



u(t,x)^Uoo, p(t,x)^p OQ as \x\ —* +oo, (1.2) 
where p^ is a given nonnegative constant. 

The term KpVD[p] corresponds to the capillarity which is supposed to model capillarity 
effects close to phase transitions in [13]. The classical Korteweg's capillarity term is 
D[p} = Ap. 

Based on Korteweg's original ideas Coquel, Rohde and their collaborators in [3] and 
Rohde in flTJ choose a nonlocal capillarity term D which penalizes rapid variations in 
the density field close from the interfaces. They introduce the following capillarity term: 

D[p] =(j)*p- p 

where 6 is chosen so that: 



^JnC^R^nW 1 ' 1 ^), f <p(x)dx = l, <f> even, and <j> > 0. 



This choice of capillarity term allows to get solution with jumps, i.e with sharp interfaces. 
Before tackling the global stability theory for the system (NSK), let us derive formally 
the uniform bounds available on (p,u). 
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1.2 Energy spaces 

We assume to simplify that P(p) = ap 1 with 7 > 1. Let II (free energy) be defined by: 

n(s)=s (r^ dz )' (l3) 

so that P(s) = sll'(s) — II(s) , n'(/>) = and if we renormalize the mass equation: 

d t U(p) + div(nn(p)) + P(p)div(u) = in P'((0, T) x R N ). 

Notice that II is convex. Multiplying the equation of momentum conservation by u and 
integrating by parts over W N , we obtain the following energy estimate: 



(1.4) 



/ (\p\u\ 2 +Il{p)+E global {p(.,t)})(x)dx{t) + f f (pD(u):D(u) 
Jrn z Jo Jr n 

+ (A + p,)\divu\ 2 )dx < / {^7T- + n(po) + E global [p Q ))dx, 

Jr n 2 P 

where we have: 

Eglobal[p(-,t)]{x) = 7 / <t>(x - y)(p{y,t) - p(x,t)) 2 dy. 

4 J R N 

The only non-standard term is the energy term E g \ oba i which comes from the product of 
u with the capillarity term KpV((f) * p — p). Indeed we have: 

k I u(t,x)p(t,x) -V([4>* p(t,-)](x) — p(t,x))dx 
Jrn 

= — k / div(u(i, x)p(t, x))([4> * p(t, -)](x) — p(t, x))dx, 
Jr n 

f 9 

= K / aiP^ * - P(t> x )) dx ' 

Jr n ot 

d f 

= ~37 / E global [p(t, ■)\{x)dx . 
at J r n 

To derive the last equality we use the relation: 
4: [ E g i oba i[p(t, -)](x)dx = 7 / [ 4>(x - y)(p(t,y) - p(t,x))-^-p(t,y)dydx 

dt J R N 2 J R N J R N at 



IR N 

where we just use integration by parts 
In the sequel we will note: 



+ 77 / / <t>(y ~ x)(p(t,x) - p(t,y))^-p(t,x)dydx, 
1 Jrn J r n at 

Jr n Jr n ot 

f d 

= -k ([4>* p(t,-)](x) - p(t,x))—p{t,x)dx . 

Jrn at 



£{p,pu)(t)= [ (lp\u\ 2 +U(p)+E global [p(.,t)])(x)dx(t), (1.5) 
Jrn 2 
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We are interested to use the above inequality energy to determine the functional space 
we must work with. 

So if we expand E g i oba i[p(.,t)](x) we get: 

E gloM [p(t,-)](x) = j(p 2 + <j)*p 2 -2p(<j>*p)). 

Because by the mass equation we obtain that p is bounded in L°°(0,T; L 1 (R N )) if we 
suppose that po G L 1 and we have supposed that 4> G L°°(R N ), we obtain that p((f> * p) 
is bounded in L°°(0, T; L 1 ^)). So we get that p 2 + <p * p 2 G L°°(0,T;L l (R N )) and as 
<fi > and p > we get a control of p in L°°(0,T;L 2 (R N )) (a property which turns out 
to be important to taking advantage of the theory of renormalized solutions, indeed p in 
L°°(0,T;L 2 (R N )) implies that p G Lf oc (R + x R N ) and we will can use the theorem of 
Diperna-Lions on renormalized solutions, see |14j ) 

In view of (11.4p . we can specify initial conditions on p/ t =o = po and pu/ t=0 = ijiq where 
we assume that: 

• po > a.e in R N , p G L 1 (R N ) n L S (R N ) with s = max(2,7), 

• tuq = a.e on po = 0, 

I 1 2 

• (defined to be on p = 0) is in L 1 ^). 
Po 

(1.6) 

We deduce the following a priori bounds which give us the energy space in which we will 
work: 

• p e L oo {0,T;L 1 (R N ) D L S (R N )), 

• p\u\ 2 e L oo (0,T;L 1 {R N )), 

• VuG L 2 ((0,T) xRT. 

We will use this uniform bound in our result of compactness. Let us emphasize at this 
point that the above a priori bounds do not provide any control on V/? in contrast with 
the case of D[p] = Ap (see in [6]). 



1.3 Notion of weak solutions 

We now explain what we mean by renormalized weak solutions, weak solutions, and 
bounded energy weak solution of problem (NSK). 

Multiplying mass equation by b (p), we obtained the so-called renormalized equation (see 

my- 

^b(p) + div(6(p)«) + (pb'(p) - b(p))divu = 0. (1.7) 

with: 

b G C°([0, +oo)) D C 1 ((0, +oo)), \b'(t)\ < ct- Xo , te{0,l], A < 1 (1.8) 
with growth conditions at infinity: 

\b'(t)\ < ct Xl , t > 1, where c > 0, -1< A a < | - 1. (1.9) 



4 



Definition 1.1 A couple (p,u) is called a renormalized weak solution of problem (NSK) 
if we have: 

• Equation of mass holds in T>'(R N ). 

• Equation holds in T>' (R N ) for any function b belonging to U.8\) and \1.9\) . 

Definition 1.2 Let the couple (po,uo) satisfy; 

• p G L l (R N ), n(p ) G L l (R N ), E global [Po] G L l (R N ), p > a.e in R N . 

• p u G {L l {R N )) N such that poK| 2 lp >o G L 1 (M Af ) 

• and such that PqUq = whenever x G {po = 0}, 

where the quantity II is defined in H^e /iaue i/ie following definitions: 

1. A couple (p,u) is called a weak solution of problem (NSK) on R if: 

(a) p G U'(U'{R N ) fors<r< +oo, 

(b) P{p) G L°°(L 1 (R N )), p > a.e in R x R N , 

(c) Vu G L 2 (L 2 (R N )) , p\u\ 2 G L°°(L 1 (IR iV )). 
fc/J Mass equation holds in V (R x R N ). 

(e) Momentum equation holds in T> (R x R N ) N . 

(f) lim^ 0+ J RJV p(t)(p = J RN p <p, V(f G V(R N ), 

(g) lim t _ 0+ f RN pu{t) ■ (f> = f RN p u ■ <p, G V(R N ) N . 

2. A couple (p, u) is called a bounded energy weak solution of problem (NSK) if in 
addition to [Id]) , J7e]) . (Tj\), [Tg] ) we have: 

• The quantity £q is finite and inequality {TT^P with £ defined by \1.5\) and with 
£q in place of £ (p(0) , pu(0)) holds a.e inR. 



1.4 Mathematical results 

We wish to prove global stability results for (NSK) with D[p] = 4> * p — p in functional 
spaces very close to energy spaces. In the non capillary case and P(p) = ap" 1 , P-L. Lions 
in [15] proved the global existence of weak solutions (p, u) to (NSK) with k = (which 
become the system of Navier-Stokes compressible isotherm) for j > ^ ^ N > 4, 7> -^rj 
if N = 2, 3 and initial data (po, uiq) such that: 



Pa, P6 



where we agree that mo = on {x G R N / po(x) = 0}. More precisely, he obtains 
the existence of global weak solutions (p, u) to (NSK) with k = such that for all 
t G (0,+oo): 
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• p G L°°(0, T; L' 1 (R )) and p € C([0, T], LP (IT)) if 1 < p < 7, 

• p G L«((0, T) x R N ) for g = 7 - 1 + % > 7. 

• p|u| 2 G L^TjL 1 ^)) and Dn G L 2 ((0,T) x E^). 

Notice that the main difficulty for proving Lions' theorem consists in exhibiting strong 
compactness properties of the density p in L p oc (R + x M. N ) spaces required to pass to the 
limit in the pressure term P{p) = ap 1 . 

Let us mention that Feireisl in [8] generalized the result to 7 > y in establishing that 
we can obtain renormalized solution without imposing that p G Lf oc (R + x R N ) (what 
needed Lions in dimension N = 2,3, that's why 7 — 1 + > 2), for this he introduces 
the concept of oscillation defect measure evaluating the loss of compactness. We refer to 
the book of Novotny and Straskraba for more details (see [IB]). 

Let us mention here that the existence of strong solution with D[p] = Ap is known 
since the work by Hattori an Li in [9], [10] in the whole space R N . In [6j, Dan chin and 
Desjardins study the well-posedness of the problem for the isothermal case with constant 
coefficients in critical Besov spaces. We recall too the results from Rohde in |17j who 
obtains the existence and uniqueness in finite time for two-dimensional initial data in 
H*(R 2 ) x H 4 (R 2 ). 

In the present chapter, we aim at showing the global stability of weak solutions in the 
energy spaces for the system (NSK). This work is composed of four parts, the first one 
concerns estimates on the density to get a gain of integrability on the density needed 
to pass to the weak limit in the term of pressure and of capillarity. The second part 
is the passage to the weak limit in the non-linear terms of the density and the velocity 
according to Lions' methods. The idea is to use renormalized solution to test the weak 
limit on convex test functions. In this part we will concentrate on the case of simple 
pressure of type P(p) = op 7 . We get the following theorem where (pm u n )neN is a suit 
of global regular approximate solution for the problem (NSK). 

Theorem 1.1 Let N > 2. Let 7 > N/2 if N > 4 and 7 > 1 else. 
Let the couple (pq,Uq) satisfy: 

• Pn is uniformly bounded in L 1 (R N ) n L S (R N ) with s = max(7, 2) and p$ > a.e in 

R N , 

• ° is uniformly bounded in L (R ), 

Pn 

• and such that PqUq = whenever x G {po = 0}. 

In addition we suppose that p° converges in L 1 (M Ar ) to pq, then up to a subsequence 
(pni u n) converges strongly to a weak solution (p,u) of the system (NSK) satisfying the 
initial condition (pq,uq) as in 111.6}) . Moreover we have the following convergence: 

• Pn^n P in C([0,T],LP(R N ))nL r ((0,T) x K) for all 1 < p < s, 1 < r < q, with 
q = s + %2 - 1 if N > 3, where K = R N except when N > 4 and 7 < f (1 + ^) 
where K is an arbitrary compact. 

• pn^nP in C([0,T],LP(R 2 ))nL r ((0,T) x K) for all 1 < p < s, 1 < r < q, with K 
an arbitrary compact inR 2 if N = 2 . 
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In addition we have: 



• PnUn -»• pu in L p (0, T; U (M )) /or all 1 < p < +00 and 1 < r < ^ 

• Pn(ui)n(uj)n ^ Pn«i«j L p (0, T; /or all 1 < p < +00, 1 < i,j < N if 
N > 3, where K = M. N except when N > 4 and 7 < + -^) where K is an 
arbitrary compact. 

• Pn(ui) n (uj) n — > p n UiUj in L p (0, T; L 1 (Q) ) /or all I < p < +00, 1 < i,j < wi/i 

an arbitrary bounded open set inM? if N = 2. 

In the third part we will focus on general pressure, and particulary van der Waal's 
pressure. In the fourth part we concentrate on the case with initial data close to a 
constant p, and we will work in Orlicz space, this case is the most adapted for the strong 
solution because it enables us to control the vacuum so that one can use the property of 
cllipticity of the momentum equation. 



2 Existence of weak solution for a isentropic pressure law 
2.1 A priori estimates on the density 

In this part we are interested by getting a gain of integrability on the density and we 
consider the case where P(p) = ap~> '. This will enable us to pass to the weak limit in the 
pressure and the Korteweg terms. It is expressed by the following theorem: 

Theorem 2.2 Let N > 2 and 7 > 1, with in addition 7 > ^ if N > 4. Let (p,u) be a 
regular solution of the system (NSK) with p > and p € L°°(L 1 nL s+£ ) where we define 
e below. Then we have if 7 > y(l + ^) for N > 4 : 

I (p 1+£ + p 2+£ )dxdt < M for any < e < —7 -1 if N > 4, 

J(0,T)xR N N 

4 

and < e < — - 1 if N = 2,3. 

with M depending only on the initial conditions and on the time T. 
7/7 < f (1 + ^) for N > 4, we have: 

[ (p 7+£ + p 2+£ )dxdt < M for any < e < ^7 - 1 */ N > 4 > 

J{Q,T)xK N 

for any arbitrary compact K with M' depending only on the initial conditions, on K and 
on the time T. 

Proof: 



We will begin with the case where N > 3 and we treat after the specific case N = 2. 
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Case N > 3: 

We apply to the momentum equation the operator (— A) _1 div in order to concentrate us 
on the pressure and we get: 

V = ^(-ArMivGH + (-A)- 1 ^.( /0 n i ^)+(2 / u + A)divu 
and in multiplying by p £ with < e < min(-^, -^7 — 1) to estimate /? 7+e , we get: 

K 



«P'- + ttP 2+£ =-^ £ (-A)- 1 div(p(V0* /O )) l/t-A)- 1 ^^,-) 



7+e 

' 2' 

+ ^(^(-AJ^divM) " [^^(-AJ^divCpu) + (fi + Qdivu 



(2.11) 

where we note £ = A + \i. We now rewrite the previous equality as follows: 

ap 7+£ + fp 2+£ = -/^(-A^divMV^p)) +p E (-A)- 1 aJ(p(« i )(« i )) 

+ — (^(-AJ^divGou)) + divtw^-AJ-Mivfcu)] + (/x + C)divu ( 2 - 12 ) 
- (p) e u • V(-A)- 1 div(pu) + (1 - e)(divu) / 9 e (-A)- 1 div(pn) , 

Next we integrate (|2.12p in time on [0, T] and in space so we obtain: 

/ (ap^ + ^-p 2+£ )dxdt = [ f|^(_A)-idiv(H] + fa + C)(div«)p £ 

+ (1 - e)(div«)p e (— A) -1 div(you) + p £ [RiRj(puiUj) - UiRiRj(puj)} 

+ div[-up £ (-A)~ 1 div(pn)] - Htp £ (- A)' 1 div (pV (<p * p)))dx dt . 



(2-13) 

where Ri is the classical Riesz transform. 

Now we want to control the term Jq J^n (p 7+£ + ^p 2+£ )dxdt. As p is positive, it will 
enable us to control II p II r 7+e and II oil r 2+ c . This may be achieved by controlling each term 

on the right side of (|2.13|) . 

We start by treating the term |j or j xRjv -^[p £ (—A)~ 1 div(pu)]. So we need to control 
p £ (-A)- 1 div(pu) in L°°(0, T; L 1 (M JV )) and pg(-A)- 1 div(p no) because: 



(0,T)xR N ^ 



^-[p £ (-A)- 1 div(pu)](t,x)dtdx = I p £ (-A)- 1 div{pu)](x)dx(t) 



Po(~ A ) div(p u )](x)dx, 



We recall that p, p 2 , p 1 and p\u\ 2 are bounded in L°°(L l ) while Du is bounded in 

2N 

L 2 ((0,T) x R N ) and it is bounded in L 2 (0, T; L~ (R N )) by Sobolev embedding. In 
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2 7 4 

particular by Holder inequalities we get that pu is bounded in L°° (0, T, (L t+i n L 5 ) (R" ) ) . 
Thus we get in using Holder inequalities and Sobolev embedding: 
p e (-A)~ 1 dw(pu) G L°°{Q,T,L l nL a ) with: 

1 e , . .7 + 1 3. 1 

— = — h mini , — ) < 1 . 

The fact that p e (— A) _1 div(pu) G L°°(0, T, L 1 ) is obtained by interpolation because 

p G L°°(L 1 ) and in using less integrability in Sobolev embedding. 

Next we have the same type of estimates for A) _1 div(/3otio)||L 1 (R JV )- 

Finally (I2.13|) is rewritten on the following form in using Green formula: 

r i ( P ^+~p 2+s )dxdt<c{i+ r I [idiH^a+K-Ar^vcpu)!] 

JO JR N 1 JO JM. N 

+ p e \RiRj(puiUj) — UiRiRj(puj)\ + K/9 e |(— A) -1 div(pV(0 * p))]dtdx). 

Now we will treat each term of the right hand side. We treat all the terms with the same 
type of estimates than P.-L. Lions in [15] , excepted the capillarity term. 

We start with the term |divu|p e |(— A) _1 div(pu)| where we have: 

|divu| G L 2 (L 2 ), p £ G L°°(L^), pu G L 2 (0, T, L r {R N )) 

with i = i + 4^r and by Sobolev embedding |(-A) -1 div(/w)| G L 2 (L S ') with ± = 
(this is possible only if r < N). We are in a critical case for the Sobolev embedding ( i.e 
r > N) only when N = 3 and 7 > 6, that's why for N = 3 and 7 > 6. 
So by Holder inequalities we get |divu|p E |(— A) _1 div(pu)| G L 1 (L Sl ) with: — = -r + | + 

i = l- I + l^<laswe have s > f . 

Moreover by interpolation |divu|p e | (— A) _1 div(pu)| belongs to L 1 (0, T; L 1 (R Ar )). 

We now treat the case N = 3 and 7 > 6 where we choose in this case e = -^7 — 1 to 

explicit precisely this case: 

|||divu| / 9 e |(-A)" 1 div(pu)||| i i < || J Du|| L 2 (i 2 ) || / 9||^ 7+e \\pu\\ 2 ( 7+E ) 6( 7+E ) 

L 7-s (L ) 
7+3 2(2 7 -3 

< 107-6 3(107-6) < C||p||^ 7+e ||pu|| 57 " 3 67 HpuHr^T^ 

L 7+3 (i 137+3 ) L2( L ^Rr) u ^ ; 

57 2(27-3) 

< r\\n\\ £ llm/ll 57-3 11/17/11 57_3) 

— IIPIlL^+ e IIP"II 67 ||P"|| 27 

L 2 (L^TE) L°°(LWT) 

^ c||pIIlt+e 

since we have i + -f- + = 1 l + + - i = 1, and ^±f) = 3 10t_6 < 3 

2 7+e 2(7+e) ' 2 7+e 6(7+e) 3 ' 57— £ 137+3 

We now want to treat the term: p 6 \(— A) _1 div(pV(0 * p))\, so we have: pV((j) * p) = 
p(V(j) * p) G L^iL 1 n L2) by Holder inequalities and the fact that we have p G L°°(L 1 ) 
and V(p G L 1 . 

After we get that /j e (-A)- 1 div(pV(0*^)) G L°°(L ri ) with: i = f + f-^ = < 1- 
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We conclude that p e {— A)~ 1 dw(pV(4> * p)) is L 00 ^ 1 ) in using interpolation when N = 
2,3. Indeed we have pV(<j) * p) G L 00 ^ 1 ) and in choosing e = j^s — 1 we have: 
1 — jj + jjS — 1 > 1 when s > + jj). This is the case when N = 2,3, and 
this is the case when N > 4 and 7 > -^(1 + jj)- 
In the other case we need to work in arbitrary compact. 

We have after the term (div(u))p e . We recall that p £ is in L°°(L~ n L~). If e > ^ (i.e 
s > \N), the bound is obvious because ^ + e > 1 and | + § < 1, we can then conclude by 
interpolation. On the other hand, this rather simple term presents a technical difficulty 
when e < 77 since we do not know in that case if divup £ E L\R N x (0,T)). One way 
to get round the difficulty is to multiply (|2.10p by P £ l{p>i}- Then we obtain an estimate 
on p s+£ l{p>i) in £ 1 ((0 > T) x R^) as p £ l {/3 > 1} |dhm| < p\divu\ £ L 1 ((0, T) x l w ) (where 
e < 5) and we can conclude since < p s+e l {p<1} < p on (0,T) x 1^ and p e L°°{L l ). 

We end with the following term p e (RiRj(p U{Uj) — UiRiRj (p Uj)) . In the same way than in 
the previous inequalities we have p £ RiRj(puiUj) is bounded in L 1 (0, T; L 1 (IR JV )). Indeed 
we have by Holder inequalities and the fact that Ri is continuous from L p in L p with 
1 < p < +00: i + 2^2 + | = i_| r + !±£<i (because s > f ). And we conclude by 
interpolation. We treat the term p e UiRiRj{puj) similarly. 

We have to treat now the case N = 2 where we have to modify the estimates when 
we are in critical cases for Sobolev embedding. 



Case N = 2: 

In the case N = 2 most of the proof given above stay exact except for the slightly more 
delicate terms p e dxvu\(— A) _1 div( / on)| and p 6 (RiRj(puiUj) — UiRiRj(puj)). 
We start with the term \p e dwu(— A)~ 1 div(pu)\. In our previous estimate it was possible 
to use Sobolev embedding on the term (— A) _1 div(pit)| only if r > N (see above the 
notation), so in the case where N = 2 we are in a critical case for the Sobolev embedding 
when 7 > 2. 

This may be overcome by using that, by virtue of Sobolev embedding, we have: 
ll|divu|p £ |(-A) _1 div(pn)||| il < C'||p||| y7+ e (Z7+E) ||p«|| 2(7+e ) 

Indeed: \ + ^- £ + ^-\ = \ + ^ < 1 = l + ^S + ^e) < 1 and l + ^+Iefej < 
1. Moreover we have as pu = ^/p^fpu 



1 



\pu\\ 2( 7 +e) < C||p||^ 7+£f/ - 7+e x 



then: 



1 



|||diH P e |(-A)-Miv(H||| L1(il) < c\\ P ^ 



L~i+^(Li+ s )- 



Next we are interested by the term p £ (RiRj(puiUj) — UiRiRj(puj)). We use the fact 
that u is bounded in L 2 (0,T;i7 1 ) and thus in L 2 (0, T; BMO). Then by the Coifman- 
Rochberg- Weiss commutator theorem in [3], we have for almost all t £ [0, T]: 
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RiRj(puiUj) - UiRiRj(puj)\\ 2(7+e) 2(7+ E ) < C\\u\\ L 2^bmo)\\p u 



2(7+e) 
L2(7+e)( L 7+e+l ) 



p £ (R i R j (pu i u j ) -UiRiRj(pUj))\\ L i < C||p|| i7 f e(L7+e) 



In view of the previous inequalities we get finally: 



/ ? IIl-t+e(Lt+ £ ) 



< C(l + ||p||? 7+s , L7+e J 




and the L 7+e (L 7+e ) bound on p is proven since ^ + e < 7 + e. 



□ 



2.2 Compactness results for compressible Navier-Stokes equations of 
Korteweg type in the case of isentropic pressure 

In the sequel we are not going to treat in details the case with N > 4 and 7 < y (1 + i), 
we just remark that the proof is the same as in the case N = 2, it suffices to localize 
because we can only apply the theorem 12.21 on the gain of integrability on any compact 



So let follow the theorem 12,21 and assume that 7 > ^- if N > 4 and 7 > 1 such that if 
(p, u) is a regular solution then p G L q ((0, T) x R^) with g = 7 + 1 — We can observe 
that in this case q > s = max(7, 2). We will see that it will be very useful in the sequel 
to justify the passage to the weak limit in some terms to get a gain of integrability on the 
density. Indeed the key point to proving the existence of weak solutions is the passage 
to the limit in the term of pressure and in the term of capillarity pX7(4> * p — p)- 
First, we assume that a sequence (p n ,u n ) nl zN of approximate weak solutions has been 
constructed by a mollifying process, which have suitable regularity to justify the formal 
estimates like the energy estimate and the theorem 12.21 (p n ,u n ) nG N has the initial data 
of the theorem 11.11 with uniform bounds, i.e: 

Moreover (p n) ^n)neN verifies the energy inequalities (|1.4j> and the previous theorem 12. 2} 
we have then: 

• p n is bounded uniformly in L°°(0, T; L 1 n L S (R N )) n C([0, T}; L P (R N )) for 1 < p < 
max(2,7), 

• Pn > a.e. p n is bounded uniformly in L q (Q,T,~M. N ) with q> s, 

• Vn n is bounded in L 2 (0, T; L 2 (R N )), p n \u n \ 2 is bounded in L°°(0, T; L 1 (IR 7V )), 

• u n is bounded in L 2 (0,T; (R N )) for N > 3. 

And we have in passing to the weak limit in the previous bound in extracting subsequence 
if necessary: 

• p n p weakly in L s ({0,T) x R N ), 



K. 
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u n ^u weakly in L 2 (0,T,H 1 (R N )), 

pl^7* weakly in L r ((0, T) x R N ) for r = 1 > 1, 

pl^>~(? weakly in L ri ((0, T) x M^) for n = | > 1. 



Notation 1 VFe wi/Z always write in the sequel B(p) to mean the weak limit of the 
sequence B(p n ) bounded in appropriate space that we will precise. 

We recall that the main difficulty will be to pass to the limit in the pressure term and 
the capillary term. The idea of the proof will be to test the convergence of the sequence 
(Pn)neN on convex functions B in order to use their properties of lower semi-continuity 
with respect to the weak topology in L 1 (M Ar ). In this goal we will use the theory of 
renormalized solutions introduced by Diperna and Lions in [?]. So we will obtain strong 
convergence of p n in appropriate spaces. 

2.3 Idea of the proof 

We here give a sketchy proof of the theorem 11.11 In this spirit we can rewrite mass 
conservation of the regular solution (p n ,Un)neN on the form: 

P ' 

— (B(p n )) + d\v(u n B(p n )) = (B(p n ) - p n B (p n ))d\\u n . 

Supposing that B(p n ) is bounded in appropriate space we can pass to the weak limit 
where we have in the energy space p n — p and u n — u, so we get: 



§- t (B(p)) + div(uB(p)) = (B(p) - P B'(p))divu (2.14) 
We will set: b{p) = B(p) - pB'{p). 

Next in seeing the mass equation for approximate solutions (p n ,u n ) n£ n , and in passing 
directly to the limit via the weak convergence arguing like P-L. Lions in [15] p 13 we get: 

^p + div(pit) = 0. (2.15) 

After we will just have to verify the passage to the limit for the product pu. Next we 
will use the theorem on the renormalized solutions of Diperna-Lions in [H] on (I2.15D in 
recalling that p G L°°(L 2 ). So we get: 

j t (B{ P )) + div(«B(p)) = 6(p)div(u). (2.16) 
Next we subtract ()2.14p to (I2.16|) . so we obtain: 

j t (Bjp) ~ B (p)) + dW(u(Bjpj - B(p))) = bjpjdbru - b(p)dwu. (2.17) 



Consequently, in order to estimate the difference B(p) — B(p) which tests the convergence 
of p n , we need to estimate the difference b(p)div(u) — 6(/?)div(u). We choose then B a 
concave function and we have: 

Bjp)-B(p)<0. 
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The goal will be now to prove the reverse inequality in order to justify that B{p n ) tends 
to B(p) a.e. 

So now we aim at estimating the difference b(p)divu — b(p)divu. This may be achieved 
by introducing the effective viscous pressure P e // = P ~ (2/U + A)divu after D. Hoff in 
[?], which satisfies some important properties of weak convergence. 

In fact owing to the capillarity term we adapt Hoff 's concept to our equation in setting: 

P eff = P + ^p 2 - (2/i + A)divn. 

Proof of theorem 11.11 

We begin with the case N > 3, and next we will complete the proof by the case N = 2 
in specifying the changes to bring. 

Before getting into the heart of the proof, we first recall that we obtain easily the con- 
vergence in distribution sense of p n u n to pu and p n (u n )i(u n )j to puiUj. We refer to the 
classical result by Lions (see |15j) or the book of Novotny and Straskraba [16]. 

Case N > 3 



We have seen in the idea of the proof [2731 that our goal is to compare B(p) and B(p) for 
certain concave functions B . From the mass equation we have obtained: 



d t (B(p) - B(p)) + div(u(B(p) - B(p))) = 6(p)div(u) - b(p)div(u). (2.18) 



So before comparing B(p) and B(p), we have to investigate the expression b(p)div(u) — 
b(p)div(u). By virtue of theorem 12.21 which gives a gain of integrability we can take the 
function B(x) = x £ , as we control for e small enough p s+£ . Our goal now is to exhibit 
the effective pressure P e ff, and to multiply it by p £ to extract divub(p) . We will see 
in the sequel how to compare it with b(p)div(u). So we focus on the convergence of the 
pressure and capillarity terms. 



Control of the term divub(p) 

So we take the div on the momentum equation satisfied by the regular solution. We get: 
—dw(p n u n ) + dfjipnU^ui) - CAdivu n + A(apl) =«div(p n (V0 * p n )) - -A(p^), 

(2.19) 

with £ = A + 2p. Applying the operator (—A) 1 to (|2,19p . we obtain: 

— (-A)~ div(p n u n ) + (-A) _1 5|(p n «) + [Cdivu n - apl - -pl\ 

= k(- A)' 1 div (p n (V4> * p n )) . 
After we multiply (|2.2U|) by p n with e that we choose enough small with e S (0, 1): 

[(p + ()divu n - apl - ^p 2 n ]p £ n = K^(-A) _1 div(p n (V(/) * p n )) 

8 ... ( 2 - 21 ) 

-Pn^(-A) 1 div ( Pn u n ) - p £ n (- A) 1 d ij (p n u l n ui). 
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So if we rewrite (|2.2ip . we have: 

[0 + C)divn„ - apl - ^p 2 n ]p £ n = «p^(-A)- 1 div(p ri (V0 * Pn )) - p e n {- A)" 1 (p^nu^) 

- Q- t ((p n ) £ {-A)- 1 div(p n u n )) + [— (p n ) e ](-A)- 1 div( /3n n n ), 

Next we have: 

[0 + C)divu„ - apl - ^pl]p £ n = rep* (-A) -1 div(p n (V0 * p n )) - p £ n (-A)~ 1 df j (p n u l n ui) 

d 

- ^[(yOn) £ (-A) _1 div(p n n n )] - div[u n (p n ) £ (- A)^ 1 div (p n u n )} 

+ (Pn) £ u n ■ V(-A) _1 div(/>„ii n ) + (1 - e)(divu n )( / 9 n ) £ (-A) _1 div(p n n n ), 

(2.22) 

or finally: 

[(ji + £) di ™n - apl - ^p 2 n ]p £ n = K P £ n (- A) _1 div(p n (V0 * p n )) 
d 

~ ^K(- A )~ div(p n u n )] - div[n n (p n ) £ (-A)- 1 div(/) n n n )] 

+ (/>n) £ [ U „.V(-A)- 1 div( /0nUn ) - (-Ar'd^pului)] 

+ (1 - e)(divu n )(p n ) £ (-A)~ 1 div(p n u n ) . 

Now like in Lions [?] we want to pass to the limit in the distribution sense in (|2.23|) in 
order to estimate divu(p) e . 

Passage to the weak limit in (12.231) 

In this goal we use the following lemma of P-L Lions in [15] to express the weak limit of 
the non-linear terms. 

Lemma 1 Let^l aopenofM. N . Let (g n , h n ) converge weakly to (g,h) in L P1 (0, T, L P2 (fi))x 
L qi (0,T, L q2 {£l)) where 1 < px,P2, Qi, Q2 < +°o satisfy, 

1 1 1 - 1 1 1 - 1 
Pi Qi P2 Q2 

We assume in addition that: 

dg n 
~dt 

and that: 



is bounded in L 1 (0,T, W m,1 (f2)) for some m > independent of n. (2.24) 



\\h n - h n (-,- + 0\\Ln(o : T,L"2(n)) as |£| -> 0, uniformly in n. (2.25) 
Then, g n h n converges to gh (in the sense of distribution on Q x (0, T)). 



14 



So we use the above lemma to pass to the weak limit in the four following non-linear 
terms of ([27Z5]) : 

T n = u n p e n (-A)- 1 div(p n u n ), Tl = p e n (-A)- 1 div(p n u n ), 

T n 3 = ^(-A)- 1 div( / 9„(V0 * Pn )), T n 4 = (divn„)(p„,) £ (-A)- 1 div(p n n ri ) . 
So we choose the different g n and h l n as follows: 



for 


- 1 n 


gl 


= U n (p n ) £ 


9 1 


= up £ 


hi = 


(-A)" 


-1 div(p n -u n ) 


for 


rp'2 

1 n 


in 


= Pn 


9 2 


= 7 


hl = 


(-A)" 


_1 div(/9 n u n ) 


for 


rp3 

1 II 


ai 


= Pn 


9" 


= 7 


hl = 


(-A)" 


_1 div(p n (V0 * p n )) 


for 


rpA 

1 n 


4 
9n 


= {d\VU n ){p n ) e 


9 4 


= divu p £ 


< = 


(-A)" 


_1 div(p n u n ). 



To show that u n {p n ) e converges in distribution sense to up 6 we apply easily lemma [T] 
with h n = u n and g n = p n . We now want to examine each term and apply the above 
lemma to pass to the limit in the weak sense. 

We start with the first term T^. We have that p e n u n G n l?(U) with ± = j- s + \ 

and ^ = ^2W^ + § = I — F + f- ^ n addition the hypothesis (I2.24D is immediately verified 
(use the momentum equation). 

We now want to verify the hypothesis (|2.25j) . so we have h\ belongs to L°°(W^? (R N )) n 
L 2 (W^'(R N )) with ^ = \ + i and A = ^1 + I = 1 + I_£. This result enables 
us to verify the hypothesis (|2.25p by Sobolev embedding. 

So we can choose (with the notation of the above lemma) q\ = 2 and q<i G (r x ' 1 



' N—r' >■> 

pi = 2, and pi = 1 — p which is possible by interpolation. Indeed we have: \ + - = 
l-A + l±,<i. 

We proceed in the same way for T 2 and T*. 

We can similarly examine Tl, because p £ n G L°°(L~nLf ) and/j n (V<£*p„) G L^^nia), 
we can choose p2 = ^, we have then (A) -1 divp n (V</> * /? n ) G L°°(0, T; '2) so that we 

N- 

can choose gi = 2, c/2 G (1, rrjg ) . We can conclude by interpolation. 

JV 2 

Finally we have to study the last non linear following term that we treat similarly as 
P-L. Lions in [15] : 

A n = ( / 9 n ) £ [-u n .V(-A)- 1 div( / 9 n u n ) - (-A)- 1 9|(p n (« i ) n (u i ) n )]. 

We can express this term A n as follows: 

A n = (p n ) e [u n ,Rij\{p n u n ). 

where Rij = (— A)~ 1 dfj with Ri the classical Riesz transform. 

Next, we use a result by Coifman, Meyer on this type of commutator (see [5]) to take 
advantage of the regularity of [u J n , Rij](p n u n )- 

Theorem 2.3 The following map is continuous for any N > 2: 

W 1,Tl (R N ) N x L T2 (R N ) — > W 1 ' V3 (R N ) N . 

(2.26) 

(a, b) -» [aj,i?ii?j]6j 



15 



To pass to the weak limit in A n we will use the previous lemma. We start with the case 
with s > 3. This quantity belongs to the space L l {W 1,q ) provided that Du n G L 2 (L 2 ) and 
pv? € L 2 (L r ) where 1 = ^ + 1 = 1- ^ + 1 in which case i = i + i= l- -^ + i<l. 

for After we can use the above lemma applied to h n = [Rij , Un](p n u l n ) and g n = p £ n . 
We can show easily in using again lemma Q] that h n converges in distribution sense to 
[Rij,Uj](p n Ui). 

So we can take: q\ = 1, p\ = +oo and Q2 G (q, jj^), P2 = 1 — ^> this one because we 

can use interpolation and we can localize as we want limit in distribution sense. 

In the case where s < 3, a simple interpolation argument can be used to accommodate 

the general case. It suffices to fix L r2 (R N ) in the application (|2.26|) and use a result of 

Riesz-Thorin. 

Finally according to the equation (|2.2ip . and after passing to the limit we get: 

[Qx + £)dh^7 - (ap~<+ £ ) - = A) -1 div(p(V0 * p)) - A) _1 div(pu)] 

- div[/^ii(-A) _1 div(pn)] +^[it.V(-A) -1 div(pu) - (-A) -1 % (pit^)] 

+ (1 - e)divu p e (-A) _1 div(/9u). 

(2.27) 



Inequality between the terms p e divu and divu p £ 

Now we are interested in estimating the term p £ divu in order to describe the quantity 
p s divu — divu p e before considering the quantity p e divu — divu p £ . We pass to the weak 
limit directly in (|2.2(jp and we get in using again the lemma ??: 



Q 

— (-A) MivO) + (-A) l d^{ P UiUj)+[{n + £)dwu - ap 



(-A)- 1 div( /O (V0* /O ) + ^p 2 . 



(2.28) 



Now we just multiply (|2.28|) with p £ and we can see that each term has a distribution 
sense. 

So we get in proceeding in the same way as before: 
[(/x + Odivup^ - {ap"i) ~{f) - = p ? (-A)- 1 div( / 9(V0* /3 )) 

-^^[pi-A^divipu)) +^[u.V(-A)- 1 div(pu) - (-A)- 1 %(p W )] ( 2 - 29 ) 



- div[p e M(-A) _1 div(pu)] + (1 - e)divu(/?) £ (-A) _1 div(pit). 
In subtracting (|2.29[) from (|2.27p . we get: 



(p + £) divu(p) e — ap"< +£ — ~^P 2+£ = (P + C)divn p e — apT /3 e — ~^P 2 P £ a - e • 
Next we observe that by convexity: 



(p7+£)7+ £ > (p£), (p7+£)7+ £ > (p7) a .e 



16 



So we get: 

divu (p) e > divu~ff. (2.30) 
i 

Comparison between p and p e e 

As on since (p n , u n ) are regular solutions we get the equality (|2.16p applied to B{x) = x £ . 
So we get: 

d 

-Q t {pnf + dw(u n (p n ) £ ) = (1 - e)dwu n (p n ) £ . (2.31) 
And after we pass to the weak limit in (|2.3ip and we get: 





—p £ + div(u p £ ) = (1 - e)divup £ . (2.32) 
at 



When we pass to the weak limit in combining with (|2.30p we get: 





— {p) £ + dW(u(p) £ ) > (1 - e)dWu(p) £ . (2.33) 



Now we wish to conclude about the pointwise convergence of p n in proving that (p £ ) e = p 
and to finish we will use the following theorem (see [8] p 34) applied to B{x) =x? which 
is convex. 

Theorem 2.4 Let (f n )neN be a sequence of functions bounded in L 1 (M Ar ) such that: 

weakly in L l (Jt N ). 

Let f : E — > (— oo,+oo] be a upper semi- continuous strictly concave function such that 
(f(v n ) G L 1 (M Ar ) for any n, and: 



(f(v n ) — ^ cp(v) weakly in L 1 (M JV ). 

Then: 

v n {y) a.e. 
extracting a subsequence as the case may be. 

Now we want to use a type of Diperna-Lions theorem on inequality (I2.33p . Our goal is to 

renormalize this inequality with the function B(x) = x^ so that one can compare p and 
i 

p £ E . Although (|2.33p doesn't correspond exactly to the mass equation, we can use the 
same technics to renormalize the solution provided that p £ L°°(L 2 ) which is the case. 
In our case it is very important that p £ L°°(L 2 ), indeed it avoids to have supplementary 
conditions on the index 7 like for the compressible Navier-Stokes system in [15] . We 
recall of Diperna-Lions theorem on renormalized solution for the mass equation. 

Theorem 2.5 Suppose that p £ L°°(L 2 ) and (3 G C^QO, oo);K). 
We have then: 

^ + divG0(p) u) = ((3(p) - p[}\p))dwu 

in distribution sense. 
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We now want to adapt this theorem for our equation (??) with (3(x) = , so we may 
regularize by oj a (with u> a = -^w(-) where u> G Cq°(M> n ), supp u> G B\ and J udx = 1) 
and find for all (3 G C§°([0, +00)): 

d — — — 

— (p £ * uj a ) + div[u p £ * Lo a ] > (1 — e)dwu p £ * oo a + R a 

where we have: 



R a = div[u p £ * uj a ] — div(u p £ ) * tu a + (1 — e) [divu p e ] * u> a — (1 — e)divn p £ * u a 
— {[3(jf * u a )) + div[n 0Qf * u) a )\ > (1 - e)dwu'p T * u a 0' (ff * uj a ) 

+ (divn) [f3(jf * u a ) - ~ff * u) a f3 Qf * u a )] + R a /3 (ff * u a ) 

= -e{dwu)lpj(3'(p^) + (divn) (3(f) + RjQf * u a ) . 

After we pass to the limit when a — > and we see that R a tends to in using lemma 
on regularization in |14j p 43. This looks like a rather innocent manipulation but it's at 
this point that we require to control p in L 2 (0, T; M. N ). And in our case we don't need to 
impose 7 > y for N = 2, 3. Hence: 

i-MW)) + div[u(3(W)} > ~ e(divuWl3'(W) + (divu){3(f). 



We then choose (3 = (^> M )~ where * M = M^(jj), M > 1, ^ G Cfi°([0, +00)), *(x) = x 
on [0, 1], supp^ C [0,2], and we obtain: 

d 1 

— (* M (^)-) + div[u¥ M G^)«] 



> (divu)^ A/ (^)t- 1 ^ M (^)^+ (divu)tf M (^) = 

> divn^M^^'^M^)- - *'m(FW\1q?>m) 

> -C |divn|Mh ( ^ >M) . 



where C = sup{|*(x)| * 1 \W(x) - xtf (x)\, x G [0, +00)}. 
Now we claim that: 



-((p^)-)+div(n(^)-)>0. (2.34) 



For proving that, we notice that by convexity {p) ££ < p, so we get : 

\\\divu\Mhj^ >M \\ L i T{L1{RN) < HdivnlliK^^llpl^^illi^^^jv) -» as M — > +00. 

We have concluded by dominated convergence. We have concluded by dominated con- 
vergence. 
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At this stage we subtract the mass equation to (12.34)) and we get in setting r = p — [p) e E '■ 

d 

— (r) + div(ur) < 0. (2.35) 

We now want to integrate and to use the fact that r > to get that r = a. a. To justify 
the integration we test our inequality against a cut-off function of the form <p(-g) where 
ip 6 Cg°(R N ), ip = 1 on 5(0,1), Suppc^ C B(0,2) and R > 1. We test the equation 
(I2.35P against ip R and we get: 

/ ^-[r(t,x)}ip R (x)-u(t,x)r(t,x)^Vip{^-)dtdx <0. (2.36) 

J[0,T]xR N Ot K K 

Next we notice that: 



f I x 

/ u(t, x)r(t, x)— V(p(— )dt dx\ < \\u\\ 2n \\r\\ 

J[0,T]xR* R R Li(0,T;LiV-2(RiV)) 11 " V (0,T;L N+2 



2JV 



* ^II Vi ^IIl°°(r^)- 



It implies that: 



f I x 

/ u(t,x)r(t,x)—Vip(—)dtdx^O as R — ► +cxd. 

J[0,T]xR n R R 



We have then: 

d 



r(t,x)ip R (x)dtdx = / r(T,x)(p R (T,x)dx - / r(0,x)ip R (0,x)dx. 

[O.TJxR^ at JR N JR N 

Now we want to verify that r(0, •) = to conclude because. Indeed we will obtain that: 
lim / r(T,x)ip R (T,x)dx — > / r(T, x)dx < and r > 0. 

then r = 0. 

We know that p n is uniformly bounded in L oc (L 1 r\L s (M. N )), then p e n is relatively compact 
in C([0,T]; LP — w) with 1 < p < s (where LP — w denote the space IP endowed with 
weak topology). Moreover (po)n converges to p^, we deduce then r(0) = a.a. 
Now as r = we conclude in using the theorem 12.41 p n converges a.a to p and that 
p n converges to p in LP([0,T] x B R ) for all p G [l,q) and in L Pl (0,T, L p ' 2 (B R )) for all 
Pi G [1, +oo), P2 G [1, s) and for all R G (0, +oo). 



Conclusion 

We wish now conclude and get the convergence of our theorem in the total space. 

We aim at proving here the convergence of p n in C([0, T], LP(R N )) n L q ' (R N x (0,T)) 

for all 1 < p < s, 1 < q < q. We have just to show the convergence of p n to p in 

C([0, T], L 1 (M 7V )). To this end, we introduce d n = ^fp~n~ which clearly converges to ^fp in 

L 2 ^{Q,T,L 2 ^{B r ))C\L 2 p{B r x (0,T)) to d = Jp for all R G (0,+oo). 

We next remark that p G C([0, T], L X {R^)) and thus d G C([0, T], L 2 ^). Indeed, 
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using once more the regularization lemma in [14) we obtain the existence of a bounded 
p a G C([0,T],L 1 (R N )) smooth in x for all t satisfying: 

+ dw(up a ) = r a in L 1 ((0, T) x R^) as a -> 0+. 

with r a = div(tt/9 a ) — div(pu) * «j a (where u; is defined as in the previous part). 
p a -» p in L 1 (R iV x (0,T)), p a / t=0 -> pA =0 in L 1 ^) as a -> 0+. 
From these facts, it is straightforward to deduce that: 

- Pv\ + div(u\p a - p v \) =r a -r rj 

and thus: 

rT 



sup / \p a — prjldx = / |r a — r^cfe. 
[0,T] Jr n ' Jo Jr n 

Since p G C([0, T], L p {Br) — w) (for all -R G (0, +00), 1 < p < s), we may then deduce 
that p a converge to p in C([0, T], L 1 (R N )). 

Next, we observe that we can justify as we did above that d n and d satisfy: 

^P + dw(u n d n ) = i(i„div(u n ), 
— — \-dw(ud) = —ddvv(u). 



Therefore once more, d n converges to d in C([0,T],L (R ) — it?). 

Thus in order to conclude, we just have to show that whenever t n G [0, X 1 ], t n — > t, 
then d n (t n ) — ► d(t) in L 2 (M 7V ) or equivalently that J RN d n (t n ) 2 dx = f RN p n (t n )dx — > n 
J RN d(t) 2 dx = f RN p(t)dx. This is the case since we deduce from the mass equation, 
integrating this equation over R N and justifying the integration exactly like previously 
that: 

p n (t n )dx = I (po)ndx I p dx = I p(t)dx. 



We then conclude by uniform continuity that ||p n (*n) — Pnif)^L l tends to 0. 
Case N = 2 

First of all, the main difficulty is the fact that we no longer have global L p bounds on 
u n . That's why most of the proof is in fact local and we know that u n is bounded in 
L 2 (0, T; L p (B R )) for all p G [1, +00), R G (0, +00). 
As we need to localize the argument, we get the following limit: 

{{p + i)diwu n -apl-^pl)p e n ^ n (( M + £)divu-a^-^)^ in V (R N x [0, T]) 



let ip G Cg°(R ), < ip < 1, supp(/3 C K for an arbitrary compact set K G R . We 
apply the operator (A) _1 div to the momentum equation that we have localized and we 
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pass directly to the weak limit: 

q 

— (-A) -1 div(<ppu) + Rij(<ppuiUj)+[(fj, + £)divu</3 - apiip] 



= k(- A)^ 1 div (ip p(V(j) * p) --ipp 2 )) + (-Ay'R. 

(2.37) 

with: 

R = dnpdj(puiUj) + (dijtp)puiUj — (p + £)A<^divu — (2p + £)Vy> • Vdivu + pAu ■ Vip 
+ Aipaf/y + aVip ■ V/P 7 + - Aipp 1 + -Vtp ■ Vp 1 . 

Now we multiply (|2.37p with p e and we verify that each term has a sense. 
So we get in proceeding in the same way as before, we can verify that (p n ) £ (—A)~ 1 R n 
converges in distribution sense to p £ (—A)~ 1 R for small enough e. We get as in the 
previous case for N > 3: 



K- 



<p[(fi + £) dwu(p) £ — ap"/ +e — — p 2+£ ] = + £)(divu)ip p £ — atppi p £ — —<pp 2 p £ ] a.e . 



We then deduce the following inequalities as in the previous proof: 
d :Q?)> + div(u(p g )«) > in £>'((0,T) x R^) 



eft 

We see that the only point left to check is the justification of the integration over M 2 
i 

of terms like div(p £e u) or div(pu) and more precisely that the integral vanishes. This 
is in fact straightforward provided we use the bounds on p G L°°(L 1 (]R Ar )) and p\u\ 2 G 
L OQ (L 1 (R JV )) and so pu G L 00 (L 1 (R JV )). Then, letting tp G C^°(M 2 ), < y> < 1, y> = 1 on 
B(0, 1) and </? = on C B(0, 2). We set </?#(■) = for R > 1, we have similarly as in 

the previous case: 

1/ eft / ru • Vtp R (x)dx\ < HV^Hloo^) — ||pit|| x i (0 T)xM 2 } -» as R -> +oo 
io JR 2 ■« 

We can then conclude as in the previous proof. 
□ 



Proof of the convergence assertion on p n u n 

We now want to show the convergence of p n u n to have informations on strong convergence 
of u n modulo the vacuum. We recall in this part some classical inequalities to get the 
convergence of p n u n , for more details see Lions in [15] We use once more a mollifier 
k a = ^k(^) where k G Cq°°(M. n ) and we let g a = g * k a for an arbitrary function g. 
We first observe that we have for all y < p < s: 

!((,*)„ -,*)W| = | / Wt,»)- fc ttx)Kft»)iW.-»)* 

+ p n (t,x)((u n ) a - U n )(t,x)\ 
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We have in using Holder inequalities with the measure k a (x — y)dy: 

f I p p-i 

\{{PnUn)a ~ PnU n ){x)\ < [ / | p n (t, y) - p n (t, x) \ P k a (x - y)dy] p (\Un\P- 1 ) p 

JR d 

+ Pn\(u n ) a - U n \(t,x). 

Hence for all t > 



W P_1 )Jlrl +\\Pn\\Lp\\(u n ) a -U n \\ p-i 

I CX J-, £ p 



< [ SUp \\p n (- + Z) - p n \\Lp} \\u n \\ 2=1 + \\p n \\Lp\\{Un)a - U n \\ p-l . 
\z\<a 'LP LP 

Next if we choose p > -Jxo, zrr < tto then ||(u n ) Q — ii ra || P -i converge to as a 

iv+z pi iv z L 2 (0,T;L~P~) 

goes to 0+ uniformly in n. In addition, the convergence on p n assure that sup| 2 | <a ||/5n( - + 
z ) ~ Pu\\lp converge to as a goes to 0+ uniformly in n. Therefore in conclusion, 
(PnU n )a — PnU n converge to in L 2 (0, T; L 1 ) as a goes to 0+ uniformly in n. 
Next (p 

n u n)a is smooth in x, uniformly in n and in t G [0,? 1 ]. Therefore, remarking 
that -§j:(pnUn)a is bounded in a L 2 (0,T; H m ) for any m > 0, we deduce that (p n u n )a 
converge to (pu) a as n goes to +oo in L 1 ((0, T) x M w ) for each a. Then using the bound 
on p n it n in L°°(L s + 1 ), we deduce that /3„n ra converge to pu in L ((0,T) x R ) and we 
can conclude by interpolation. 

The last convergence is a consequence of the strong convergence of p n and p n u n . 



3 Existence of weak solution with general pressure 

In the sequel we concentrate us only on the cases iV = 2, 3. We now want to extend 
our previous result to more general and physical pressure laws. In particular we are 
now interested by two cases, the first one concerns monotonous pressure law (close in a 
certain sense that we will precise to p 1 pressure) , the second one is the case of a slightly 
modified Van der Waals pressure. 

The technics of proof will be very similar to the previous proof, only technical points 
change. 

3.1 Monotonous pressure 

In this section, we shall investigate an extension of the preceding results to the case of a 
general monotonous pressure P, i.e P is assumed to be a C 1 non-decreasing function on 
[0, +oo) vanishing at 0. 

We want here to mention in the general situation our new energy inequality, we recall 
the inequality (|1.4|) : 

I (lp\u\ 2 + U(p)+E global [p{. 1 t)])(x)dx(t)+ [ I {pD{u):D{u) 
JR N 1 JO JR N 

+ (\ + p)\dwu\ 2 )dx < [ (B- + IL(p ) + E global [p ])dx. 

Jrn 2p 
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where we define II by: ^ ) = -sr for all t > 0. 

There are two cases worth considering: first of all if P(t) is such that Jq ^§^ds < +00 
then we can choose H(p) = p f£ ^r-ds. 

In the other case, i.e lim^o = c > 0, we can choose H(p) = p ^r-ds and IT behaves 
like plog p as p goes to 0. 

We now consider a sequence of solutions (p n , u n ) and we make the same assumptions on 
this sequence as in the previous section except that we need to modify the assumptions 
on p n . We assume always that p n is bounded in C([0,T],L 1 (R )), P{p n ) is bounded in 
L oo {0,T,L 1 (R N )). 

■ (Pn)n>i is bounded in L°°(0, T, L 2 ), 
and we also assume that we have: 

(p £ n P(p n ))n>l is bounded in L\K x (0,T)) 
for some e > 0, where K is an arbitrary compact set included in M> N . 

Theorem 3.6 Let the assumptions of theorem \l.l\ be satisfied with in addition P a mono- 
tone pressure. 

Then there exists a renormalized finite energy weak solution to problem (NSK) in the 
sense of definitions li.il and li.M Moreover P(p n ) converges to P(p) in L l (K x (0, T)) 
for any compact set K . 

Proof: 

In this situation, the proof of theorem 13.61 is based on the same compactness argument 
as in the theorem ll.il In particular, there is essentially one observation which allows us 
to adapt the proof of theorem 11.11 Namely we stills obtain the following identity for the 
effective viscous flux: 

f3(.Pn)((p + 6divn n - P(p n ) - ~pl) ^ n J(p){{p + £)divu - F\p) - |p 2 ), 
with (3{p) = p £ for e small enough. We have then: 

(p + - P{p)p e - 1 7^) = (P + £)div^ - Pjpp - ^ff), (3.38) 

Now we can recall a lemma coming from P.-L. Lions in [15] : 

Lemma 2 Letpi, P2 S C([0, 00)) be non- decreasing functions. We assume thatpi{p n ), 
P2(Pn) o-re relatively weakly compact in L l (K x (0, T)) for any compact set K C M. N . 
Then, we have: 

Pi(p)P2(p) > Pi(p)P2(p) a.e. 

We get finally as in the proof of theorem 11.11 in using lemma [2] 

divu p £ > divnp e , 

All remaining argumentation of the proof of theorem 11.11 can be performed to conclude. 

□ 
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3.2 Pressure of Van der Waals type 

In this section we are interested by pressure of type Van der Waals which consequently 
are not necessarily non-decreasing. That's why in the following proof we will proceed 
slightly differently. 
So we consider the pressure law: 

RT p 

Pip) = -. — ap 2 for p < b — 9 for some small 9 > 

b-p 

and we extend the function P to be strictly increasing function on p > b — 9. 
We have then: 

1. P' is bounded from below, that is: 

P'(p) > ~P for all p > 0. 

2. P is a strictly increasing function for p large enough. 

Under the above conditions, it is easy to see that the pressure can be written as: 

P{p)=P 1 {p)-P 2 {p). 
with Pi a non-decreasing function of p, and 

P 2 G C 2 [0, +oo), P 2 > 0, P 2 = for p > p. 



Remark 1 The a priori energy estimate give us the bound of p in L°°(L 2 ), then we have: 

T\\p\\ 2 

\{(t,x) G (0,T) x R N /\p(t,x)\ > b}\ < ^f 2(R }) . 

Then we can controll in measure the set where P is different from the Van der Waals 
measure, and it's a set of finite measure. 

We obtain the following theorem. 

Theorem 3.7 If in addition to the above assumptions, we assume that p^ converges in 
LMjH N ) to po then (p, u) is a weak solution of the system (NSK) satisfying the initial 
condition and we have: 

p n p in Ci[0,T], L p (R N ) f] L r i(0,T) x R N ) for all 1 < p < 2, 1 < r < 1 + — , 



Proof: 

Most of the proof of theorem 13.71 is similar as theorem 11.11 We will use a approxi- 
mated sequel T^ (introduced by Feireisl in [?]) of p by some concave bounded function. 
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Definition 3.3 We define the function T E C 00 ^) as follows: 

T(z) = z for z€ [0,1], 

T(z) concave on [0, +00), 

T(z) = 2 for z > 3, 

T(z) = -T{-z) for z£ (-oo,0], 



AndT k is the cut-off function: 



T k (z) = kT{- k> 



In following the proof of theorem ll.lt we get: 

d 

g^i L k{p) ~ L k {p)) + dw((T k (p) - T k {p))u) + T k {p)dwu - T k (p)dwu 



= {T k (p) - T k (p))divu in V ((0,T) x 

where L k (p) = plog(p) for < p < k, and < L k (p) < plog(p) otherwise. 
So we get in integrating in time on [ti,^ 



(L k (p)-L k (p))(t 2 )dx- / (L fc (p) - L k (p))(tx)d. 

rt 2 r rt2 



+ / / T k (p)divu — T k (p)divu dxdt = / / (T k (p) — T k (p))divudxdt. 
Jti Jr n Jti 

We can show that: 

\\ T k(p) ~ Tk(p)\\L 2 ({0,T)xR N ) ~^k- 



For poving the previous inequality, we see that \\T k (p) — T k (p) Hl^^tJxR^) — > f° r 
k — > +00. We then conclude by interpolation with T k (p) — T k {p) € L q (0,T) x W N ) with 
q > 2. By Holder inequality we obtain that: 

/ / (T k {p) — T k (p) )divu dxdt — > for A; — > +00. 

Jt! Jr n 



We have then: 



lim / (L k (p) - L k (p))(t 2 )dx - [L k {p) - L k {p))(ti)dx 

k^+QO J R N J R N 

rt 2 



k 



lim / / T k {p)d\\u — T k {p)di\u dxdt. 



(3.39) 



We set: 



dft[ Pn ^ p](t) = lim / (L k (p)-L k (p))(t)dx 

k^+OO J R N 

A(k,p)= / / T k {p)d\vu — T k {p)divu dxdt. 
Jh Jr n 
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We can show as in the previous proof of theorem 11.11 that : 

ft2 r 

/ T k {p)d\vu — T k {p)d\vu dxdt 
ti J K 

*2 



lim I" f ((P(p n ) + ^pl)T k (p n )-(P(p) + ~p2)T^pj)dxdt. 

n^+oo J h J K 2 2 



for any compact K C ~S, N . 

Using the lemma [2] we deduce that: 

rt 2 



hm / 2 / (PM + ^p 2 n )T k {p n ) - (P 1 (p) + %f?)T k (p)dxdt) < 0. 

n-*+ocJ tl J RN 2 2 

We have then: 

dft[p n -► p](t 2 ) - dft[p n -> p](ti) < lim ( lim f [ P 2 {p n )T k (p n ) 

fc-»+oo \n-++oo J tl J RN 



-P 2 {p)T k (p)dxdtj. 

As the sequence (/3 n )neN is bounded in L°°(0, T, L 2 (M. )), and P2 is a bounded function, 
we have: 

lim f lim f [ P 2 (p n )T k (p n ) - P 2 {p)T k (p)dxdt] 

\n^+oo J Ti J RN J 

= lim / \ P 2 (p n )p n - P 2 (p) pdxdt. 
J. J RN 



Since the function P 2 is twice continuously differentiable and compactly supported in 
[0, +00), there exists A > big enough such that both p — > Aplogp — pP 2 (p) and 
/> — > Aplogp + pP 2 (p) are convex functions of p, indeed the twice derivative are positive. 
As a consequence of weak lower semi-continuity of convex functionals, we obtain: 



lim f f P 2 (p n )p n - P 2 (p) pdxdt 
J. J RN 



<A / (plogp- plogp)dxdt+ / (P 2 (p) - P 2 {p)) pdxdt 



Futhermore we have: 

,f 2 



{P 2 {p)-P 2 {p))pdxdt< / / (P 2 (p)-P 2 (p))p 



< A / / (p logp — p log p)dxdt 

Jtl J p<p r 

< Ap r {p l°g P — P l°g p)dxdt 

Ju Jr n 



The previous relation gives: 



ft2 

dft[ Pn p](t 2 ) < d/t[p„ p](t) + W / d/t[p n -» /£>](*)• 

ti 
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Applying Gronwall's lemma we infer: 

dft[p n -> p](t 2 ) < dft[p n -> p](ti)exp(w(t 2 - ti)). 
We conclude that dft[p n — > p](t) = Vi, because p° converges strongly in L 1 to po- 

4 Weak solutions with data close to a stable equilibrium 

We consider in this section one situation which is rather different from the three cases 
considered in the preceeding sections. This situation is relevant for practical applications 
and realistic flow and they involve conditions at infinity different from those studied. 
We wish to investigate the system (NSK) with hypothesis close from these of the theorem 
for strong solutions. We want then to study the system with a density close from a stable 
equilibrium in the goal to can choose initial data which avoid the vacuum. We look now 
for a solution (p, u) defined on R x R N of the system (NSK) (where P(p) = ap 7 ) with 
p > on R x R N . 

In addition we require (p, u) to satisfy the following limit conditions: 
(p,u)(x,t) — > (p, 0) as \x\ — ► +00, for all t > 

where p > 0. 

Such an analysis requires the use of the Orlicz spaces. We define the Orlicz space Lp(R N ) 
as follows: 

L«(R N ) = {/ G L} oc (R N )/fl {lfm e L P (R N ), /l {m > 5} G L«(R N )}. 
For the following properties, we can give this definition too. 

Definition 4.4 We define ^ as a convex function on [0, +oo) which is equal (or equiv- 
alent) to x p for x small and to x q for x large. 

L«(R N ) = {/ G L} oc (R N )/*(f) e L\R N )}. 

We can check that Lp(R N ) is a linear space. Now we endow Lp(R N ) with a norm so that 
Lp(R N ) is a separable Banach space. 

ll/llL ? (R-)=mf(t>0/ *({)<1). 
We recall now some useful properties of the Orlicz space: 
Corollary 1 We have: 

1. Embedding: 

L q p (R N ) C L q p \(R N ) if l< qi <q< +oo, 1 < p < pi < +oo. 

2. Topology: f n -^ n in L p (R N ) if and only if^(f n ) and that: 

\\j\\L q p (R N ) 
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We recall now some useful properties of the Orlicz space: 
Proposition 4.1 We then have: 

• Dual space: If p > 1 then (L g p (R N ))' = L q \ (R N ) where q = ^ p = JB^. 

• Let F be a continuous function on R such that F(0) = 0, F is differentiable at 
and F(t)\t\-° -> a / at t -> +oo. T/ien if q>6, 

F(f)£L}(R N ) if f£Ll(R N ). 

Our goal is now to get energy estimate. We have to face a new difficulty. Indeed 
p, p\u\ 2 , p 1 need not belong to L . 

We first want to explain how it is possible to obtain natural a priori bounds which 

correspond to energy-like identities. 

Next we write the following formal identities: 

1 d -{p< -p>- jp y ~ 1 (p -p)+ div[n^-(^ - pr-ip)] = u ■ V{pT) 



7 — 1 dt ' 7 — 1 

d \u\ 2 \u\ 2 ( 4 - 40 ) 
p— — - — h pu ■ V— pA.u ■ u — £Vdivu • u + au ■ V p 1 = KpuV(cj) * p — p). 

OjV — ■ Zi 



We may then integrate in space the equality (|4.40p and we get: 
( / P l ^ + -^( P ^+(7-l)p r -lP r - 1 p)+E global \p-p}dx)(t) 

JWL N 1 7 — 1 

+ fdsf 2p\Du\ 2 + 2£|divu| 2 dx < / pj^f- + -Ar + (7 - l)f ~ lf~ l Po) 
Jo Jr n Jr n 2 7 ~ 1 

+ EgldballPO ~ P}dx. 

(4.41) 



Notation 2 In the sequel we will note: 

h(p) = P 1 + (7 - !)P 7 - 7P 7 ~ V 
Now we can recall a theorem (see [15]) on the Orlicz space concerning this quantity: 
Theorem 4.8 j 7 (p) £ L 1 (R N ) if and only if p - p G L%, 

In following this theorem and our energy estimate we get that p-p e L co (0,T;Ll(R N )) 
for all Tel. 
Moreover we have: 

EglobalW, ■) - p](x) = ^(p-p) 2 + (/)*(p-p) 2 - 2(p - p) (0 * (p - p)). 

Then in using the fact that p-pe L°°(0, T; L7,(R N )) ^ L 2 (R N ) + L~f(R N ) and inter- 
polation on V0, we get that p - p G L°°(0, T; L 2 (R N )). 
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We may now turn to our compactness result. First of all, we consider sequences of solu- 
tions (p n ,u n ) of the system corresponding to some initial conditions (/o°,it°). 
In using the above energy inequalities, we assume that j 7 (/4J), -E 9 iofea«[Pn ~~ P] an d ^nl^nl 2 
arc bounded in L°°(L 1 (IR 7V )) and that 

Pn~ P converges weakly in I? 2 (M. N ) to some po — p. 

We now assume that: 

j 7 (p n ), E g i oba i[p n - p], p n \u n - u\ 2 are bounded in L oo (0,T,L 1 (R N )), 
Moreover we have for all , T G (0, +00) and for all compact sets K C W N : 

Pn-pe L°°(L 2 (R N ) and p n is bounded in L"((0,T) x K), 
for some q > s. 

Du n is bounded in L 2 (R N x (0,T)) 

u n is bounded in L 2 (0,T, H 1 (B R )) for all R, T G (0, +00) . 

Extracting subsequences if necessary, we may assume that p n , u n converge weakly re- 
spectively in L 2 ((0,T) x B R ), L 2 {Q,T;H l {B R )) to p, u for all R,T G (0,+oo). We also 
extract subsequences for which p n u n , p n u n (g> u n converge weakly as previously. 

Remark 2 We notice that the conditions at infinity are implicitly contained in the fact 
that (p n - p) 2 and p n \u n \ 2 G L 1 (R JV ). 

We then have the following theorem. 

Theorem 4.9 Let 7 > 1. We assume that p\ converges in L X {B R ) (for all R G (0, +00) ) 
to po- Then (p n ,Un)n^N converges in distribution sense to (p,u) a solution of (N SK). 
Moreover we have for all R,T G (0, +00): 

p n ^ n p in C([0,T},L p (B R x (0,T)) D L Sl (B R x (0,T)) for all 1 < p < s, 1 < s a < q. 
with q = s + jj — 1 . 
Proof: 

As in the theorem ll.il we want test the strong convergence of p n on concav function B. 
Since the proof is purely local, we have again for small enough e > 0: 

(pnY{{p + £)divu n - a(p n y - ^p 2 n ) -^ n (p) £ ({p + £)divu - djfi - ^p 2 ) 

in V'((0,oo) x R N ), 

^(p T ) + div(up T ) > (1 - e)(divu)/F in X>'((0,oo) x R N ). 

(4.42) 

Next since (p £ )~ G L 2 (B R x (0, T)) for all R,T G (0, +00), as in the theorem II. II in using 
a result of type Diperna-Lions on renormalized solutions, we get: 

+ div(u(^p) > in T> ((0, 00) X R N ) . (4.43) 
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while we have (p) £E < p a.e in x (0,+oo) and {p) £ J t=Q = P/t=o hi R N ■ 

Now in subtracting the second equality of (14. 43ft from the mass equation and setting 
i 

/ = p — p £ E , we have: 

j t (f) + div(«/) < 0, / > a.e, f /t=0 = in R N . (4.44) 

Next we want again to show from (|4.44j) that / = 0, in integrating (|4.44|) and in using 
the fact that / < to conclude that / = 0. The difference with the proof of theorem ll.il 
is to justify the integration by parts as we work in different energy space. 
We need a cut-off function. We introduce ip G C^°(M. N ), < ip < 1, suppy? C B2, cp = 
1 on B\ and we set ip R = </?(;§) for R> 1. Multiplying (|4.44|) by <pr{x), we obtain: 

4/ fVR{x)dx= [ 1/ U -V^). (4.45) 
at J r n J r n K K 

Next, if T > is fixed, we see that supp Vy(^) C {R < \x\ < 2R} , therefore, for R 
large enough we have: 

4 / fVR.{x)dx = [ -fu.Vip(^), 
at . ton . ton R R 



< 
~ R 



( , r ( 4 - 46 ) 



/ f\ u \ l (R<\x\<2R)dx, forte (0,r), 
To conclude that / = 0, we only have to prove that: 

Si fHl(R<\x\<2R)dx^0 as # +00. (4.47) 

We now use that fact that / G L°°(0, T, L 2 (R N )) and f\u\ 2 E L°°(0, T, L 1 ^) for all 
T £ (0, +00) to control (|4.47j) . The second fact is obvious since < p and p\u\ 2 £ 
L°°(0,T; L^iR 1 *)). 

1 

In order to prove the first claim, we only have to show that (p) £e -p£ L°°(0,T, L 2 (R N )). 
We rewrite {p n f - {pf = (p+(p n - p)) £ - {pf is bounded in L°°(0, T, L 2 {R N )) in using 
proposition 14.11 with F{x) = (p + x) e — (p) £ . So we have yff € ^(L 4 ^)) and we get: 

— / dt f\u\h R <\ x \< 2R )dx < -£meas(C(Q,R,2R))4. 

ft JO JR N ft 



We recall that: 
Then we get: 



meas(C(0,R,2R)) ~b_» +00 C{n)R 



N 



d_ 

Tit 



fip R (x)dx -^R^ + oo 



and we conclude as in the proof of theorem 11.11 

At this stage, it only remains to show that, for instance, p n converge to p in C([0, T], ^{Br)) 
for all R,T £ (0, +00). In order to do so, we just have to localize the corresponding ar- 
gument in the proof of theorem 11.11 
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Therefore we choose for R,T G (0,+oo) fixed, <p G C^°(R ) such that = 1 on 
< ip on Then, we observe that we have: 

Ql& 2 Pn) + div(u n (v? 2 p n )) = p n u n ■ Vtp 2 , -q^p) + div ( u (^ 2 />)) = P u ■ 
d 1 

-^(fy/Pn) + div(u n (tpp n )) = -(dwu n )(py/p^ + yfpZu n ■ Vtp, 



qI&Vp) + div(u((pp)) = -(dwu)ip^p + ^fpu ■ V99. 

From these equations, we deduce as in the proof of the previous theorem, that ip 2 p £ 
C([0,+co),L 1 (R N )), ip^/p € C([0,+oo),L 2 (R JV )) and that converges weakly in 

L 2 (R N ), uniformly in t € [0,T]. Therefore, in order to conclude, we just have to show 
that we have: 



whenever t n G [0, T], t n — > n t, and this is straightforward since we have, in view of the 
above equation: 
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